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Influence of adatom interactions on second layer nucleation 
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We develop a theory for the inclusion of adatom interactions in second layer nucleation occurring 
in epitaxial growth. The interactions considered are due to ring barriers between pairs of adatoms 
and binding energies of unstable clusters. The theory is based on a master equation, which describes 
the time development of microscopic states that are specified by cluster configurations on top of 
an island. The transition rates are derived by scaling arguments and tested against kinetic Monte- 
Carlo simulations. As an application we reanalyze experiments to determine the step edge barrier 
for Ag/Pt(lll). 

PACS numbers: 68.55.-a, 68. 55. Ac, 81.15.-z, 81.15.Aa 



Whether thin films in epitaxy become rough or grow 
smoothly layer by layer depends on the onset of sec- 
ond layer nucleation on top of islands in the first layer: 
If the rate of this nucleation is large, it is likely 
that mounds are formed before layer completion, while 
small second layer nucleation rates favor layer by layer 
growth. Detailed theories have been developed in the 
past 0, |[ ||, ||, |j| to predict in dependence on the island 
size R, the ratio a s = (v a /v t ) exp(— AE s /k B T) of the 
hopping rates v t exp(— Ejj/k B T) and v s exp(—E s /k B T) 
on top of the island and over the island edge, respec- 
tively, and the ratio D/Fa 4 of the adatom diffusivity 
D = (v t a 2 /4) exp(—Eu/k B T) and deposition flux F; a is 
the lattice constant. Studies so far focused on adatoms 
that form stable nuclei once more than i of them cluster 
together, but otherwise do not interact. 

In this Letter we show how adatom interactions 
strongly influence the second layer nucleation. We con- 
sider two types of these interactions: (i) Clusters of size 
smaller than i+1 are metastable, i.e. the detachment of 
an adatom or the break up of clusters into two pieces re- 
quires a dissociation energy E dis > ksT, and (ii) the ap- 
proach of two atoms by diffusion is hindered by additional 
barriers. As a generic model we consider a barrier E IiDS 
in form of a ring with radius £ around each adatom. The 
occurrence of such barriers was explored by Fichthorn 
and Scheffler based on extensive density functional cal- 
culations 0| . Similar interaction effects can be caused by 
Shockley surface states Q . 

The aim of this Letter is to develop a theory based 
on rate equations that allows us to take into account 
these interaction effects for nucleation processes in con- 
fined geometry as encountered on top of a growing island. 
As an application we quantify possible errors when de- 
termining the step edge barrier E s || in second layer 
nucleation experiments under the assumption of non- 
interacting adatoms. For the Ag/Pt(lll) system we com- 
pare the results at different levels of sophistication of the 
underlying theory for determining the step edge barrier. 

In order to take into account the adatom interactions, 
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FIG. 1: States and the corresponding transition rates involved 
in the second layer nucleation with ring barrier interaction 
and a critical nucleus of size i = 1. The adatoms in configu- 
ration 2 are inside each others' ring barrier, which is marked 
by the connecting line. The two states with configuration 
label v = 3 contain a stable dimer. 



we refer to our rate equation approach for second layer 
nucleation outlined in |}|. In this approach one consid- 
ers one compact island with radius R that evolves in time 
according to some growth law R = R(t). To clearly sepa- 
rate the different types of interaction effects, we consider 
the presence of either metastable clusters or the ring bar- 
rier. 

We denote by p n ,v(t) the probability to find the island 
in a state where in total n atoms are on top of the island 
in configuration v at time t. The configuration label v 
specifies the way the n atoms are decomposed into clus- 
ters and which pairs of atoms have a distance smaller 
than £ and are thus "weakly bound" by their ring barri- 
ers. 

In Fig. |l] the most important states are shown together 
with the rates W n ^^ n > M > connecting them The prob- 
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abilities p n ,u{t) obey the master equation 

dp n „ 



(If 



^ \W^n' : u' — >n,vPn',v' l^n,^ — >n' .i>' Pn,v\ j (1) 



where W n ^ v ^ n i y is the rate from state (n, v) to (n',v') 
and can depend on time through the growth law R = R(t) 
(see below). Equations (0) are numerically solved subject 
to the initial condition p n .u — S n fi- 

The rates W n ^^ n ' y are expressed in terms of ele- 
mentary rates Wf, Wi, W t ™l, Wf, Wb and W a (cf. 
Fig. Q). The first two rates Wf and Wj refer to pro- 
cesses not involving the ring barrier: Wf = ttFR 2 is 
the deposition rate of adatoms onto the island and Wi = 
(D / R 2 )[nia/ (a s R) + H2\~ 1 is the loss rate of adatoms 
that leave the island by surmounting the step edge bar- 
rier j|. The coefficients K\,K2 are of order one ||. In or- 
der to account for the reduced mobility of pairs of atoms 
bound via the ring barrier compared to individual atoms, 
the loss rate for such pairs is approximated to be equal 
to that of a single atom. 

The latter four rates involve the ring barrier in 
form of the associated Boltzmann factor a ring = 
exp[-AE xins /k B T], AE IiDS = E I1DS - E D , and the ring 
radius £. Wf — ttF^ 2 is the deposition ("flux") rate into 
a ring. W a = H a 2D/ (7r£ 2 ) is the attachment rate for two 
atoms in a circle which are confined to a separation dis- 
tance smaller than £. It results from the time cx £ 2 /2D 
for the two atoms to encounter by diffusion on top of an 
island with radius £. For the relevant case a < £ <C R we 
can further derive the rates for formation and breakup of 
pairs weakly bound by their ring barriers. The breakup 
rate is given by 



W b = K b 



(2) 



It results from the probability cx 27r£/7r£ 2 for two atoms 
to have distance £ times the rate cx a Ilae D / a 2 to overcome 
the ring barrier. The rate for formation of an adatom pair 
weakly bound by their ring barriers is 
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if in total n adatoms are on top of the island. The associ- 
ated time (Wri^g) -1 for two atoms results from two contri- 
butions: The time cx R 2 /D for first reaching the interac- 
tion distance £ <C R and the time for overcoming the bar- 
rier. The latter is given by [(27r£a/7ri? 2 ) x (a tins D / a 2 )] -1 , 
where the term 27r£a/7ri? 2 is due to the probability for a 
pair to have a distance in the interval (£, £ + a) and the 
term a Tias D/a 2 is the rate to overcome the ring barrier. 

~ (2) 

The rate W r Y a ' s for a single pair has to be multiplied by 
the number n(n— 1)/2 of distinct pairs to obtain the total 
rate for n adatoms. 

The validity of all formulae for the elementary rates is 
tested against kinetic Monte-Carlo (KMC) simulations 
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FIG. 2: Scaling plots of the rates for formation and breakup 
of atom pairs bound by their ring barriers. The symbols in 
the upper part of the figure refer (i) to n = 2, £/a = 3 for 
R/a = 10 (+), 22 (*), 46 (o), 100 (x), (ii) ton = 2, £/a = 6 
for R/a = 46 (□), 100 (o), and (m) to n = 3, R/a = 6 for 
£/a = 3 (a) and 6 (V). The solid line in the upper part is a 
fit to f(x) = (kix + K2) -1 [cf. eq. (H)] and the solid line in 
the lower part has slope one [cf. eq. (B)]. 



performed on a hexagonal lattice for the (111) surface. 
As an example we show in Fig. scaling plots for the 
two rates W rin g and Wb demonstrating the good agree- 
ment of the KMC data with the predictions of eqs. (PP)- 
All constants k in Eqs. (|||) are of the order of one p0[. 
The approach outlined above provides a general frame- 
work to treat the problem of second layer nucleation and 
may be applied and extended to a variety of situations. 
It furthermore allows one to gain detailed insight into 
the microscopic pathways followed during the nucleation 
process. 

To apply our theory we consider the determination 
of step edge barriers in second layer nucleation exper- 
iments. In these experiments one measures the prob- 
ability f[R(t)\ that a stable cluster has nucleated on 
top of the island until time t (i.e. the fraction of "cov- 
ered islands"). In systems with a substrate mediated 
ring barrier, second layer nucleation will be aggravated 
and nucleation sets in later. A repulsive ring barrier 
therefore has an effect similar to a reduced step edge 
barrier, and thus yields an apparent measurement value 
AE^ smaller than the "true" AEs • On the other hand, 
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metastable dimers facilitate the formation of a stable nu- 
cleus, leading to values AEg larger than A£g. 

To estimate the significance of the interactions, we first 
present results for the relative error [AE$ — AE^]/AEs 
as a function of the interaction parameters. By solving 
Eq. (Q) we obtain the curves f[R(t);AE s ] that refer to 
processes including interactions. On the other hand, we 
can solve the rate equations neglecting the interaction 
(A£ ring = 0), thus obtaining f Q [R(t); AE^]. By fitting 
the curves /o to the "true " curves / we obtain the ap- 
parent AE'g '. For the ratio of attempt frequencies Vs/vt 
for adatom hopping over the step edge and on a terrace, 
we use the generic value of one (for a general discussion 
on attempt frequencies, see jll)). Using a generic growth 
law R(t) oc Vt (with the prefactor (D/Fa 4 )^ 2 ^ 2 ^ \f~Fa? 
from standard nucleation theory p2|) we show in the 
upper half of Fig. || results for the relative error as a 
function of AE Hns . In the lower half of the figure we 
include the results of an analogous analysis for i = 2 
with metastable dimers (for a setup of the correspond- 
ing rate equations, see Because the dependence on 
D/Fa 4 is only weak, this parameter is not varied in the 
plots. We see that for small AE S , even weak interactions 
AE[ nt lead to a large relative error. Closer inspection 
shows, that for fixed ratio AE[ nt / AE S , the relative er- 
ror decreases with increasing interaction energy (cf. the 
dashed line in Fig. ||). In the case of metastable dimers, 
the error reaches a plateau when the dissociation energy 
becomes so large that a dimer is stable on the time scale 
of the formation of a stable trimer. This signifies the 
transition to the i = 1 case. 

For Ag/Pt(lll) detailed second layer nucleation mea- 
surements of the type discussed above were performed by 
Bromann et al. |13[| . The system Ag/Pt(lll) is partic- 
ularly suited as a reference, since it has the advantage 
that many of the relevant parameters were determined 
both by experiment lO, O, [l5, 16 and by first-principle 
calculations |t], [Tt], [19). Experimentally, the diffu- 
sion barrier Ejj for Ag adatom diffusion on strained Ag 
islands grown on Pt(lll) is E D ^ 60 meV 61. It has 
been shown by density functional calculations that sil- 
ver adatoms diffusing on top of an already existing silver 
island on the platinum surface exhibit a strong barrier 
A-E r i n g = -E r ing — En = 50 meV at a distance £/<z = 2.1. 

The original analysis by Bromann et al. was based 
on a mean field-type theory for second layer nucleation 
("TDT approach", see 6J). However, it has been shown 
recently that even in the absence of interaction effects 
it is necessary to re-analyze the data with an extended 
theory fi || (see also Qj). In this theory fluctuation- 
dominated regimes occur for critical nuclei i < 2 ^ , and 
it yields, for the experimental conditions and a prefac- 
tor i/ t = 10 9 Hz used in 6j|, the second-layer nucleation 
rate Q(R) = ^KxFR 5 /(Ta s ) where T = D/Fa 4 . By 
contrast to the generic growth law R(t) oc \ft used for 
Fig. |^, in the analysis of Bromann et al. one has to deal 
with an exponential growth law R(t) (for details see [pi[). 
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FIG. 3: Relative error of the apparent additional step edge 
barrier AE^ when neglecting interactions (AE- in t = 0) with 
respect to the "true" barrier AE S calculated for A_Bi nt > 0; 
upper part: AEi nt = AE rills , lower part: AEi nt = Ai? d i a - The 
symbols refer to different values AEs/ksT = 3 (•), 6 (■), 9 
(a), 12 (o), and 15 (□). The dotted line in the upper figure 
is drawn for AE S = AE TiDS . All plots are calculated with 
D/Fa 4 = 10 9 and the generic growth law R(t) oc yt. 



Metastable states are not considered here, because i = 1 
in the temperature range of the experiment. In Fig. ^ we 
have fitted the measured fraction of covered islands f(R) 
yielding AE^ = 65 meV, which is about twice the value 
obtained previously based on the TDT approach (see Ta- 
ble I). For the much less reliable prefactor ratio J2(| we 
find Vs/vt = 41. The rate equation theory without inter- 
actions yields essentially the same value AE^ = 68 meV 
thereby supporting the scaling approach [^|, |[. 

The next step in theoretical sophistication is the inclu- 
sion of the ring barrier. When using the value AE Iiag — 
50 meV [Q, we find a significantly decreased value of 
AE S — 52 meV. This decrease seems to be in contradic- 
tion to the reasoning above that AE^ < AE S . However, 
the values were obtained by fitting both the prefactor 
Vs/vt and AE S to the data and cannot be directly com- 
pared because the prefactors turn out to be very different 
(see table I). When fitting under the constraint i^/^t = 1 
W, the value for AE S increases from 42 meV to 48 meV 
in accordance with the general argument. Interestingly, 
the behavior of the standard error (see table I) indicates 
that the system is better described when including the 
interactions. However, this interpretation should not be 
driven too far, since the standard error is only an indica- 
tion of the quality of the fit, but should not be interpreted 
as the measurement error. 

It was recently argued that the often reported small 
prefactors for hopping on weakly corrugated surfaces are 
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due to neglecting the interaction effects in the analysis of 
first layer nucleation experiments |Tl[ . We therefore also 
analyze the data using the theoretical prediction = 
10 12 Hz jf9| in the rate equation approach. The three a s 
values obtained for each temperature shown in Fig. ^ lie 
very closely on a line in an Arrhcnius plot (see the inset 
in Fig. ^). Both when using u s /vt as fitting paramter 
and when setting v s /i>t = 1, we find comparatively large 
values AE S — 74 meV and AE S = 82 meV, respectively. 
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FIG. 4: Fraction of covered Ag islands of size R as measured 
by Bromann et al. Jl3[ ]. The solid lines refer to fits with the 
scaling approach when neglecting interactions and the dashed 
lines refer to fits with the rate equation theory when including 
the ring barriers. The inset shows the Arrhenius fit for values 
a B obtained by the rate equation theory for us = ft = 10 12 Hz. 
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AE S [meV] 
(v s /vt = l) 


10 9 





TDT 


30 ±5 
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30 


10 9 





Scaling 


65 ±5 


41 


43 ± 2.2 
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Rate Eq. 


68 ±1 


83 


42 ±2.6 
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50 


Rate Eq. 


52 ±1 
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48 ± 0.4 


10 12 


50 


Rate Eq. 


74 ± 2 


0.3 


82 ±0.8 



TABLE I: Comparison of results for the additional step edge 
barrier for Ag/IML Ag/Pt(lll) obtained with theories of dif- 
ferent levels of sophistication (see text). The standard errors 
from the fitting procedure are shown for the values of AE S . 



In summary we have developed a rate equation theory 
that allows us to quantify the influence of certain interac- 
tion effects on the nucleation on top of islands in epitaxial 
growth. The formalism is rather general and may be ex- 
tended to other types of adatom interactions and other 
kinds of geometrical constraints. We applied the theory 
to second layer nucleation experiments designed to ex- 
tract the additional step edge barrier. It was shown that 
the value of AE S depends very sensitively on whether the 
interactions are carefully accounted for in the analysis. 

We are indebted to W. Dieterich and F. Scheffler 
for very helpful discussions and thank the Deutsche 
Forschungsgemeinschaft for financial support through 
the SFB 513 and the Heisenberg program (P.M., Ma 
1636/2-1). 
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